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Qs ^ We review the present knowledge of disk accretion in young low mass stars, 

and in particular, the mass accretion rate M and its evolution with time. 
The methods used to obtain mass accretion rates from ultraviolet excesses 
and emission lines are described, and the current best estimates of M for 
Classical T Tauri stars and for objects still surrounded by infalling envelopes 
are given. We argue that the low mass accretion rates of the latter objects 

fT^ ' require episodes of high mass accretion rate to build the bulk of the star. 

^SJ ' Similarity solutions for viscous disk evolution suggest that the inner disk 

mass accretion rates can be self-consistently understood in terms of the disk 
mass and size if the viscosity parameter a ~ 10~^. Close companion stars 
^ ' may accelerate the disk accretion process, resulting in accretion onto the 

Ijr^ ' central star in < IMyr; this may help explain the number of very young stars 

fT^ , which are not currently surrounded by accretion disks (the weak emission T 

fT^ , Tauri stars). 
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I. INTRODUCTION 

The initial angular momenta of star-forming molecular cloud cores 
i-G ' must be responsible for the ultimate production of binary (and mul- 

tiple) stellar systems and circumstellar disks. Unless the protostellar 
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Q ' cloud core is very slovifly rotating, much or most of the stellar mass is 

pi , likely to land initially on a disk, and must subsequently be accreted 

^ ■ from the disk onto the protostellar core. In the early phases of this 

accretion, the circumstellar disks may be relatively massive in compar- 
ison with the central protostar, and so gravitational instabilities may 
K> , be important in angular momentum transport and consequent disk evo- 

^ • lution. At later phases, disk evolution is likely to be driven by viscous 

5t 1 processes, perhaps limited by condensation of bodies. Stellar magnetic 

fields may ultimately halt the accretion of material onto the central 
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star. There are substantial uncertainties in our theoretical understand- 
ing of these processes, and we must rely on observations for guidance. 
In this chapter we review the present knowledge of the disk accre- 
tion process around low mass stars, in particular the rate at which mass 
is transferred onto the star, M, and its evolution with time. Knowl- 
edge of M will help put lower limits on the disk mass at a given age, 
independent of uncertainties in dust opacities. It puts constraints on 
disk physics, and in particular on temperatures and surface densities 
at a given age, and thus on conditions which obtain during the time 
when solid bodies agglomerate. We begin with a description of the 
different methods of determining disk accretion rates and the values of 
M obtained applying these methods to stars in different groups and 
environments. We then discuss the implications of the observational 
results for disk accretion physics and evolution, and in particular in- 
vestigate whether the data are consistent with simple models of viscous 
disk evolution. 



II. DETERMINATION OF M 

A. Infrared excesses 

When the infrared excess emission in Classical T Tauri stars (here- 
after CTTS) was recognized as emission from dust at low tempera- 
tures distributed in a circumstellar disk, it was thought that the ex- 
cess energy would be a direct measurement of the accretion luminosity 
{Lace = GMM^/R^,, where M* and i?* are the stellar mass and ra- 
dius). However, it soon became apparent that their spectral energy 
distributions (SEDs) did not follow the law XF\ oc X"^^^ expected for 
standard accretion disks, but were much flatter (Rydgren and Zak 1987; 
Kenyon and Hartmann 1987). It has now been realized that the major 
agent heating the disks of many CTTS is irradiation by the central star 
(Adams and Shu 1986; Kenyon and Hartmann 1987; Calvet et al. 1991, 
1992; Chiang and Goldreich 1997; D'Alessio et al. 1998). The optically- 
thick disks of CTTS are probably "flared", i.e., have "photospheres" 
that curve away from the disk midplane; this makes them more effi- 
cient in absorbing light from the central star; this extra heating helps 
to increase the disk scale height and thus increases the flaring. Self- 
consistent calculations (in various approximations) indicate that the 
flaring is especially important at large radii (Kenyon and Hartmann 
1987; Chiang and Goldreich 1997; DAlessio et al. 1998). This makes 
it extremely difficult, if not impossible, to say anything about accretion 
energy release, and thus accretion rates, in outer disk emission. 

For many CTTS it is very difficult to extract quantitative estimates 
of mass accretion rates even from emission from the inner, flat parts of 
the disk. The basic reason is that the accretion luminosities of many 
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CTTS are smaller than the luminosity the optically-thick inner disk 
produces as a result of absorbing light from the central star. 

The effective temperature T of the disk is determined by internal 
viscous dissipation and external irradiation by the central star; at a 
given annulus, 

T\r)^Tj + T,\ (1) 

where T^ = [{3GM^,M /Snar^)/] is the effective temperature that 
would result from accretion without irradiation, with / = 1 — ^JJlL^Jr)^ 
and Ti^ = (2/37r)T,'*(_R*/r)'^ is the effective temperature for the case of 
only irradiation (the flat disk approximation). In these expressions, we 
assume that the mass accretion rate through the inner disk is constant 
and equal to the rate at which mass is transferred onto the star. The 
mass accretion rate at which T^ ~ T„ is 

calculated for typical Taurus CTTS parameters (K7-M0, age ~ 1 Myr) 
and r ^ 3i?*. Since the median value of M in CTTS is - 10'^ MQyr~^ 
(see below) , this implies that irradiation dominates the inner disk heat- 
ing for a large number of the stars and determines the amount of flux ex- 
cess. Only for objects with significantly high M will the (near-infrared) 
disk emission depend upon the mass accretion rate. 

Other considerations affect a disk's near-infrared emission. It was 
originally thought that CTTS disks extended all the way into the star, 
and disk material joined the star through a "hot boundary layer" , where 
half of the accretion energy was dissipated, while the other half was 
emitted in the disk (cf. Bertout 1989). It has now become apparent 
that for typical values of the stellar magnetic field (< 1 kG, Basri et 
al. 1992) and of the mass accretion rate, the inner disk of CTTS can 
be disrupted by the magnetic field (Konigl 1991; Najita et al., this vol- 
ume); material falls onto the star along the magnetic field lines, forming 
a magnetosphere, and merges with the star through an accretion shock 
at the stellar surface. In support of this model, fluxes and line profiles 
of the broad (permitted) emission lines in CTTS are well explained by 
the magnetospheric flow (Calvet and Hartmann 1992; Hartmann ct al. 
1994; Edwards et al. 1994; Muzerolle et al. 1998a,b,c; Najita et al., 
this volume), while the ultraviolet and optical excess fluxes are well 
accounted for by the accretion shock emission (Calvet and Gullbring 
1998). 

For truncation radii Ri ~ 3 — 5 i?, , the maximum temperature in 
a disk with typical parameters will be ~ 1200 — 850K, so that disk 
emission drops sharply below r^ 2 — 3.5/im. The calculation of the disk 
luminosity from the observed excess also depends upon the cosine of 
the inclination angle, which is generally unknown. As an illustration of 
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these points, Meyer et al. (1997) have shown that the so-called "CTTS 
locus" in JHK and HKL diagrams, i.e., the region populated by CTTS 
outside that corresponding to reddened main sequence stars (Meyer et 
al., this volume) could be explained in terms of emission by irradiated 
accretion disks with inner holes of different sizes (but inside co-rotation 
radii), and a random distribution of inclinations. However, no one-to- 
one correlation could be found between the excess and M. In a sample 
of stars in Taurus with known reddening and M, the largest excesses in 
the locus were produced by stars with the highest M, but the converse 
was not true, due to the effects of holes and inclination. 

Finally, in objects surrounded by infalling dusty envelopes, emis- 
sion from the dust destruction radius, which peaks at ^ 2/ini, may con- 
tribute significantly to the K band and longward (Calvet et al. 1997), 
hiding the intrinsic disk emission. 

The difficulty of deriving reliable measures of M from excess in- 
frared emission requires using alternative methods of estimating M's. 
We will consider such methods in the next two sections. 

B. Veiling luminosities 

The best evidence for CTTS disk accretion generally comes from 
the interpretation of the ultraviolet and optical spectra. The photo- 
spheric absorption lines are "veiled" , that is, they are less deep than 
standard stars of the same spectral type. This veiling is produced by 
(mostly) continuum emission from a region hotter than the stellar pho- 
tosphere. The luminosity of the veiling or excess continuum in CTTS 
is typically ~ 10 % of the total stellar luminosity. It is difficult to ac- 
count for this extra energy with a stellar origin, and impossible in the 
case of most extreme CTTS, where the excess continuum luminosity 
is several times the stellar luminosity. This conclusion is reinforced by 
the existence of the weak-emission T Tauri stars or WTTS. The WTTS 
have similar masses and ages as the CTTS in many regions, but do not 
exhibit the UV-optical veiling continuum of CTTS, showing that the 
excess emission is not an intrinsic property of young stars. The veiling 
emission is observed only when there is excess near-infrared emission 
(Hartigan et al. 1995; Najita et al., this volume), strongly supporting 
the idea that accretion from a disk is occurring, producing both the 
infrared excess as well as the hot continuum, as originally envisaged by 
Lynden-Bell and Pringle (1974). 

In the magnetospheric model, the disk is truncated in the inner 
regions, which is precisely where an accretion disk would emit most of 
its energy. The total disk accretion luminosity is expected to be about 

Lacc(disk) - GMj'I/{2Ri) + Ld,ss 

where Ldiss is any energy that might be dissipated in the disk by the 
stellar magnetic field lines (Kenyon et al. 1996). If we assume that the 
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effect of the stellar magnetic field is to substantially reduce the angular 
momentum of the disk material, so that it starts out nearly at rest at 
Ri before falling onto the star along the magnetic field lines, then the 
infalling material should dissipate its energy at the stellar surface at a 
rate ~ {GM^,M /R^){1 — R^,/Ri). For disk truncation radii ~ 3 — 5i?*, 
most of the accretion energy is released at the stellar surface in the hot 
accretion shock whose radiation is observed as the veiling continuum 
(Konigl 1991; Calvet and GuUbring 1998). 

Accretion luminosities can be obtained from measurements of the 
veiling of the absorption lines by the following procedure in which excess 
and intrinsic photospheric emission components are separated. The ob- 
served fluxes are fltted by the scaled, dereddened fluxes of a standard 
star, assumed to be of the same spectral type as the object star, plus a 
continuum flux, which produces the observed veiling at each absorption 
line wavelength. This flt yields the spectrum of the excess continuum 
and the reddening towards the star. The total excess luminosity is cal- 
culated from the measured luminosity using a model to extrapolate the 
emission to unobserved wavelengths. Finally, the mass accretion rate 
can be calculated from stellar mass and radius estimated from the lo- 
cation of the star in the HR diagram and comparison with evolutionary 
tracks. 

Accretion luminosities for CTTS in Taurus have been determined 
from veiling measurements by a number of authors, including Harti- 
gan et al. (1991), Valenti et al. (1993), and Hartigan et al. (1995). 
The mass accretion rates derived from these studies range from ~ 
lO"^M0j/r~^ for Hartigan et al., to ^ lQ^^MQyr~^ for Valenti et al. 
More recently, GuUbring et al. (1998, GHBC), re-measured accretion 
rates for a sample of Taurus CTTS using spectrophotometry covering 
the blue and Balmer jump region of the spectrum down to the atmo- 
spheric cutoff, and found their results to agree with the previous lower 
estimates. In a sample of 17 CTTS in Taurus, GHBC found a median 
value of M ~ lO~^M0j/r~^, with a factor of ~ 3 in estimated error. 

These differences reflect the cumulative effects of a number of dif- 
ferences in working assumptions rather than a fundamental difference 
in approach. Among them we can list the following: (1) Adopted evo- 
lutionary tracks; (2) Adopted model of accretion. The magnetosphere 
model predicts more luminosity per mass accretion rate than the old 
boundary layer model; (3) Adopted physics. Some treatments assumed 
that a substantial fraction of the emitted accretion energy was absorbed 
by the star, which seems unlikely (Hartmann et al. 1997); (4) Differing 
methods used to determine extinction corrections. The photospheres 
of T Tauri stars show color anomalies relative to standards, which ren- 
der the determination of reddening uncertain. In particular, GHBC 
found large color anomalies in some non-accreting WTTS often used 
as templates; this would cause the extinction to be overestimated in 



6 CALVET ET AL. 

all objects. These anomalies are probably due to spots on the stellar 
surface and/or unresolved, cooler companion stars. 

To provide a powerful tool to determine mass accretion rates for 
large samples of stars for which short wavelength spectrophotometry 
is difficult to obtain, GHBC determined the relationship between the 
accretion luminosity and the excess luminosity in the U photometric 
band, Ljj, for the stars in their sample. Figure la shows the excellent 
relationship between log Lace and log Lu- A least-squares fit to the 
line yields: 

log (Lacc/Le) = 1.04^^8 log (Lu/Lq) + 0.98t^:l]? (2) 

The spectrophotometric sample from which equation (2) was derived 
is formed by stars in the K5-M2 range, with most of the stars be- 
ing K7-M0. The application of this calibration to a wider spectral 
type/mass range remains to be confirmed, but theoretical models of 
accretion shock emission from stars of differing mass and radii indicate 
that for the characteristic energy fluxes found in the accretion columns 
of CTTS, the spectrum of the excess emission does not depend on the 
underlying star, and the proportion of the total excess luminosity in the 
U band (~ 10%) displayed by equation (2) holds (Calvet and GuUbring 
1998). 

Mass accretion rates have been determined for a larger sample of 
TTS in Taurus using equation (2) by Hartmann et al. (1998, HCGD). 
The median mass accretion rate is ^ 10^^ MQyr~^, similar to that of 
the spectrophotometric sample, but the errors are larger, given the lack 
of simultaneity in the photometric measurements and the high degree 
of variability of CTTS. HCGD also used spectral types and photometry 
from the compilation of Gauvin and Strom (1992) to estimate mass ac- 
cretion rates for K5-M3 stars in the Chamaeleon I association. The me- 
dian mass accretion rate in Chamaeleon I is '-^ 4 x IQ^^MQ-yr^^. A his- 
togram of the mass accretion rate for TTS in Taurus and Chamaeleon 
I determined from the ultraviolet excess is presented in Fig. 2a. 

C. Magnetospheric emission lines 

The measurement of mass accretion rates from the optical and ultra- 
violet excess fluxes, either spectrophotometrically or from broad-band 
photometry, is sensitive to reddening corrections. In heavily-extincted 
stars, the UV fluxes are either very uncertain or unobservable; thus, in 
early stage of evolution such as the infall/protostellar phase, or in very 
young, dense regions of star formation, such methods cannot be used. 
For these objects, it is necessary to devise methods to measure M in 
spectral regions that are less affected by intervening dust. 

Most emission lines present in the spectra of young objects are 
thought to be produced in the magnetospheric flow (Najita et al., this 
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volume). Since the material flows through the magnetosphere with a 
rate similar to M in the inner disk, the emission fluxes of the lines 
formed in this flow are expected to depend upon the mass accretion 
rate. Theoretical models show this to be the case, but other parameters, 
such as the unknown temperature structure and the characteristic size 
play a role too, as well as the optical depth of the line (Muzerolle 
et al. 1998a). For these reasons, empirical correlations between line 
luminosities and accretion luminosities have been investigated, leaving 
the interpretation to future theoretical work. 

Muzerolle et al. (1998b, c, MHCb,c) have undertaken spectroscopic 
studies in the red and infrared of the sample of stars with accretion rates 
determined from spectrophotometric measurements, and have found 
remarkably good correlations between the luminosity of the Ca II triplet 
8542, Pa/3, and Br7 lines with accretion luminosities. Figure lb show 
the correlation between Lace and the luminosity in Br7, as an example. 
Least-squares fits to the data give: 

log {Lacc/Lrs) = (0.85 ± 0.12) log (Lca/78542/i0) + (2.46 ± 0.46) , (3) 

log [Lacc/Lo) = (1.03 ± 0.16) log (Lpap/LQ) + (2.80 ± 0.58) , (4) 
and 

log (Lacc/Le) = (1.20 ± 0.21) log (^s.^/i©) + (4.16 ± 0.86) . (5) 

These subsidiary calibrators of the accretion rate provide the means to 
determine accretion rates in heavily extincted objects for the first time. 

MHCc used luminosities of the Br7 line for extincted Class II 
sources in p Oph from Greene and Lada (1996) to make the first de- 
termination of accretion luminosities for these sources. They find that 
the distribution of accretion luminosities is similar to that in Taurus. 
Using approximate spectral types from Greene and Meyer (1995), the 
estimated median mass accretion rate is ~ 1.5 x 10^^MQyr~-^, although 
the span of spectral types covered is wider than in the Taurus sample. 

Even more far reaching, MHCc determined accretion luminosities 
for the deeply embedded Class I objects (sources still surrounded by 
infalling envelopes) for the first time, using Br7 luminosities. Since 
typically Ay ^ 20 — 30 for Class I objects, significant extinction is 
expected at Br7 {Ak ^2 — 3), which introduces large uncertainties in 
the determination of the line luminosities. The light from the central 
star-|-disk and from the envelope itself is absorbed and scattered by 
the infalling envelope, in proportions that depend on uncertain param- 
eters such as the geometry of the inner envelope (Calvet et al. 1997). 
Following Kenyon et al. (1993b), MHCc calculate a correction factor 
K-KO, by assuming that the central object has J-K colors similar to 
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CTTS, and estimating the intrinsic J magnitude from the bolonietric 
luminosity. 

Figure 3 shows accretion luminosities estimated by MHCc for a 
small number of Class I sources in Taurus and in p Oph which have data 
on Br7 (Greene and Lada 1996) plotted against the total luminosities. 
The accretion luminosities of Class I sources are significantly lower 
than the system bolometric luminosity; in fact, the mean accretion 
luminosity is ~ 10 - 20% of the mean bolometric luminosity of the 
sample. This fact implies that the luminosity of Class I is dominated by 
the stellar component, and gives a natural explanation to the similarity 
of the distributions of luminosities of Class I and optically visible T 
Tauri stars in Taurus (Kenyon et al. 1990). 

Determining the mass accretion rate from Lace requires knowledge 
of M^:/R^,. We can obtain this ratio assuming that the stars are still on 
the birthline (Stabler 1988), which seems justified since they are still 
accumulating mass from the envelope/disk. We used the observed bolo- 
metric luminosity to locate the object along the birthline (calculated 
for an infall rate of 2 x lO~^M0yr"^). Figure 2b shows the distribution 
of M thus obtained for Class I sources in Taurus and p Oph. 

The mean mass infall rate in the envelope, M^, has been estimated 
in Taurus from fitting the spectral energy distributions and scattered 
light images of Class I sources to be ~ 4 x lQ~^MQyr~^ (Kenyon et al. 
1993a), a value which is consistent with theoretical expectations (Shu 
et al. 1987). The present determination of disk accretion rate shown 
in Fig. 2b indicates that many Class I objects are slowly accreting 
from their disks, despite the fact that mass is being deposited in the 
outer disk at a higher rate. This discrepancy was first recognized by 
Kenyon et al. (1990), as the so-called luminosity problem of Class I 
sources. If infall were spherical, the luminosity of the system would be 
given by an accretion luminosity ^ GM^,Mi/R^,. But in this case, the 
predicted accretion luminosity should be of order IOLq, while the mean 
luminosity of Taurus Class I sources is ~ 1 — 2Lq (Kenyon et al. 1990, 
1994). Kenyon et al. (1990) pointed out that the mass accretion rate 
in the disk onto the star does not necessarily equal the mass infall rate 
from the envelope onto the disk, since they are regulated by different 
physical processes. The imbalance between the infall and accretion 
rates in this picture leads to an accumulation of mass in the disk. 
These disks could eventually become gravitationally unstable (Larson 
1984), with consequent rapid accretion until sufficient mass has been 
emptied out of the disk. Kenyon ct al. (1990) suggested that these 
episodes could be related to the FU Orionis outbursts. 

D. The FU Ori disks 

The FU Orionis outbursts are now recognized as a transient phase 
of high mass accretion in the disk around a forming low mass star 
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(see review by Hartmann and Kenyon 1996 and references therein). 
Although so far FU Ori objects have been mostly studied as isolated 
phenomena, it is becoming increasingly clear that episodes of high mass 
accretion rate may be a crucial, if not dominant, process in the forma- 
tion of stars. We briefly review here the determination of Ad for these 
objects. 

The canonical FU Ori objects were discovered from their increase 
in brightness by several magnitudes over time scales of months to years 
(Her big 1977), during which the luminosity increased from values typ- 
ical to CTTS to a few hundred Lq. Since the emission is dominated 
by the accretion disk in FU Ori objects, the accretion luminosity can 
be readily determined from the observed SED. Additional information 
is required to independently obtain M^^j R^, and M, which can be es- 
timated from the surface temperature in the inner disk (i? ~ 2 - 3 
R*)-, Tmax ~ 7000K{Lacc/W0LQf^'^{R^/2RQf^^ (from T„ in equation 
[1]), and measurements of the rotational velocity. Using this method, 
M ~ lO~^M02/r~^ has been inferred for the canonical objects for which 
outbursts have been observed, consistent with their high luminosities, 
> IOOLq. 

A significant number of objects have been identified as FU Ori 
objects in recent years, mostly at IR wavelengths. Photometric vari- 
ability indicative of outbursts has been detected in only a few ob- 
jects; their classification as FU Ori systems has been based mostly 
on the presence of the near-infrared first overtone bands of CO in 
deep absorption, comparable to late type giants and supcrgiants and 
to the canonical FU Ori objects. Most of these objects are in very 
early stages of evolution, as is their prototype L1551 IRS 5, embed- 
ded in infalling envelopes and associated with Herbig-Haro objects 
jets and/or molecular outflows (Elias 1978, Graham and Frogel 1985 
Reipurth 1985; Staude and Ncckcl 1991; Kenyon et al. 1993; Han- 
son and Conti 1995; Hodapp et al. 1996; Reipurth and Aspin 1997 
Sandell and Aspin 1998). The luminosities of the objects identified 
so far as members of the FU Ori class range from ^ IOLq to 8OOL0 
Since M - 2 x 10-'^MQyr-'^{Lacc/'i-0LQ)x[{M^/R^)/0.18MQ/RQ]~^ 
this range implies M > few xlO~^ MQyr~^ (assuming they are on the 
birthline, with typical M^,/R^, ~ 0.18Mq/Rq). This value is consistent 
with the presence of CO in absorption, since disk atmosphere calcula- 
tions indicate that the temperature of the continuum forming region is 
higher than the surface temperature (neglecting any wind contribution) 
for M > IQ-'^MQyr-^ (Calvet et al. 1992). 

The nature of the low luminosity FU Ori objects and their relation- 
ship to the canonical objects remains to be elucidated. One possibility 
is that disks undergo instabilities driving outbursts of different magni- 
tude. Alternatively, these objects could be in the phase of decay from a 
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canonical FU Ori outburst. An argument in favor of the second possi- 
bility comes from comparison of the mass loss rate required to drive the 
molecular outflow and the present mass accretion rate. For L1551 IRS5 
and PP13S {L ~ IOLq and 30L©, respectively), the momentum flux of 
the molecular flow is ~ few xlO~''M02/r~^kms~^. (Moriarty-Schieven 
and SncU 1988; Sandell and Aspin 1997). For typical velocities of the 
jet of few xlOOkms"^, and assuming momentum conservation, the mo- 
mentum flux implies mass loss rates of the same order to the inferred 
mass accretion rates, while both theoretically and observationally, the 
ratio between the mass loss rate and mass accretion rate is close to 
~ 0.1 (Calvet 1997). This discrepancy may imply that the mass accre- 
tion rate of the disk was much higher when the material giving rise to 
the molecular outflow was ejected. 

III. THE EVOLUTION OF ACCRETION 

A. Observed M vs. age 

Figure 4 shows mass accretion rate vs. age for CTTS in Taurus, p 
Oph, and Chamaeleon I, for which significant infall from an envelope 
has stopped. We also show the range of M covered by Class I sources 
in Taurus, assuming a median age of 0.1 Myr (estimated from the 
ratio of the number of Class I sources to T Tauri stars in Taurus, and 
adopting a mean age of 1 Myr for the latter, Kenyon et al. 1990). The 
data indicates a clear trend of accretion decaying with time (HCGD; 
see also Hartigan et al. 1995), even in the relatively short age spread 
covered by the observations. There is also a very large spread in mass 
accretion rate at a given age, which makes it very difficult to quantify 
the rate of decay. HCGD shows that the slope of a least-squares-fit to 
the data is highly dependent on the errors of M and age. If errors are 
assumed larger in M than in age (the most likely situation, cf. HCGD), 
then the slope is w —1.5 (with large uncertainty). 

B. Disk Masses 

We can use the observed mean decay of M with time to obtain an 
estimate for the mass of the disk for comparison with masses estimated 
from dust emission. The amount of mass accreted by a CTTS from the 
present time to infinity constitutes a lower limit to the mass remaining 
in the disk. If M{t) ~ M{to){t/toY'^ , then this limit to the disk mass 
is Mace = M{to)to/{v - 1) -- 2M{to)to, with rj -- 1.5. This gas mass 
can be compared to the disk mass estimated from dust emission in the 
submillimeter and millimeter range, Mj^rm which is very dependent on 
the assumed opacities (Beckwith et al. 1990; Osterloh and Beckwith 
1995, masses corrected by a factor of 2.5, HCGD). The comparison (for 
single stars) yields log{2M x age/Mmm) = — 0.07±0.21, indicating that 
masses inferred from the current mass accretion rates and age estimates, 
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which are probing the gas, are consistent with disk masses estimated 
independently from mm- wave dust opacities, which in turn suggests 
that the dust opacities used in the latter estimates are appropriate. 

C. Early stages 

Figure 4 indicates that the mean mass accreted onto the star dur- 
ing the CTTS phase is ^ IO^^Mq. This suggests that by the time 
stars reach the optically visible CCTS stage, the remaining disk mass 
is relatively low and little mass is added to the forming star. This pic- 
ture implies that the bulk of stellar accretion must occur during the 
(highly-extincted) infall phase, when the disk is being continually re- 
plenished by the collapsing envelope. However, during the infall phase, 
our results suggest that disk accretion rates during quiescent phases are 
only slightly higher than those in the CTTS stage, precluding addition 
of more than ~ O.OIM© during these phases. Significant additions to 
the mass of the growing protostar must therefore come via material 
accreted during transient episodes of high accretion. 

Andre et al. (1993) suggest that the Class objects are in the 
earliest phases of envelope collapse, when central accretion rates are 
expected to be the highest, and argue that the Class sources are the 
true protostars. These very heavily extincted objects have somewhat 
higher luminosities than the mean Class I luminosity in Taurus, and so 
may have higher accretion rates onto the central star. However, Class 
sources tend to be less frequent than Class I sources, especially in 
Taurus; thus, the lifetime of the Class phase may be too short to 
account for most of the stellar accretion. 

The high accretion rate episodes required to explain the accretion 
of the stellar mass can be attributed to instabilities in quiescent low 
M disks, which result in outbursts and transient periods of high mass 
accretion rates. Several models have been presented to trigger those 
outbursts. The most accepted model is that of thermal instabilities in 
the inner disk (Lin and Papaloizou 1985; Clarke et al. 1990; Kawa- 
zoe and Mineshige 1993; Bell and Lin 1994; Beh et al. 1995). This 
model can naturally explain the occurrence of outbursts during the in- 
fall phase because it requires a high background mass accretion rate 
from the outer disk, of the order of a few xlO~^MQyr^^ , to match 
the observations. Moreover, according to this model outbursts will be 
triggered in the disk as long as mass is deposited in the outer disk at 
this rate, ensuring that mass will be transferred from the initial cloud 
into the star. 

The number of currently-known FU Ori objects is insufRcient to 
explain the formation of typical stars mostly through outbursts (Hart- 
mann and Kenyon 1996). However, the known sample of FU Ori disks 
is very incomplete, because the identifying characteristic of high mass 
accretion rate disks is the near-infrared bands of CO in absorption and 
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many objects may be too heavily extincted to be detected at ^ 2.2/im. 
More and more sensitive observations of embedded objects in the near- 
infrared are necessary to test this hypothesis. 

D. Viscous evolution in TTS disks 

Once the main infaU phase is over and the envelope is no longer 
feeding mass and angular momentum into the disk, we expect disk 
evolution to be driven mostly by viscous processes, namely, those in 
which the angular momentum transport is provided by a turbulent 
viscosity. In this and the next section, we attempt to interpret the 
observed properties of CTTS disks in terms of viscous evolution, with 
the aim of understanding the main physical processes at play and the 
role of initial and boundary conditions. Similarly, physical models for 
disk evolution help us relate properties and evolution of the inner disk, 
as measured by M{t), to those of the outer disk, such as radius and 
mass, Rd and Md- 

The disk angular momentum is 

Jd ^ dUnR^ oc MdR]/^. (6) 

Jo 

where Vl is the Keplerian angular velocity. If we neglect the (small 
amount) of angular momentum being added to the star, or the possible 
angular momentum loss to an inner disk wind (Shu et al. 1994; Shu 
et al., this volume), then the disk angular momentum is approximately 
constant. In turn, this requires Rd oc Md~ , so angular momentum 
conservation implies that the disk expands as the mass of the disk 
is accreted to the star. Evolution occurs on the viscous time scale, 
tvisc ~ R^/i^, where v is the viscosity (Pringle 1981). If i^ oc R'' , then 
dRd/dt ~ Rd/t^^sc oc Rd''-\ so Rd 'oc t^/(^~^\ Md oc ^-1/2(2-7)^ and 
M oc i~(5/2-7)/(2-7)^ Therefore and in principle, from the observed 
decay of M with time, M{t) oc t~^, we can obtain 7 = [27] — 5/2)/ {rj — 
1), and the evolution of radius and mass of the disk can be predicted. 
As Lynden-Bell and Pringle (1974) showed, analytic similarity so- 
lutions describing the evolution of disk properties exist for the case 
of power-law viscosity (see also Lin and Bodenheimer 1982). These 
analytical similarity solutions have been applied by HCGD as a first 
approximation to the evolution of T Tauri disks. (More complex models 
have been used by Stepinski (1998) to consider similar issues; whether 
the observational constraints justify approaches with more assumed 
parameters is not clear.) HCGD argue that the use of a power-law 
viscosity can be justified on approximate grounds. Using the a pre- 
scription for the viscosity (Shakura and Sunyaev 1973), ly = aCgH oc 
Cs^ /^{R) oc T{R)R^/^^ where Cs is the sound speed and H the scale 
height. If r(i?) oc R^"^, then v oc R^^'^^''. With q ^ 1/2, corresponding 
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to irradiated disks at large distances from the star (Kenyon and Hart- 
mann 1987; D'Alessio ct al. 1998), and also found by empirical fitting 
to apply to most disks in CTTS (Bcckwith ct al. 1990), then 7 ~^ 1, 
which is roughly consistent with the observed slope of the M vs. age 
data, 7] « 1.5 (section III. A). 

HCGD calculated similarity solutions for viscous evolution for a 
range of initial conditions applicable to CCTS disks. Figures 5a and 
5b show the evolution of mass accretion rate and disk mass for a sub- 
set of models in this study. Initial disk masses have been taken as 
Md{0) = O.IM©, consistent with small disk masses remaining after 
disk-draining episodes of high M (III.C). Values for other values of 
Md{0) can be obtained by simple scaling. Model results are shown for 
three values of the initial radius, 1, 10, and 100 AU, which cover an 
order of magnitude in angular momentum, consistent with the spread 
of angular momentum between half of the binaries in the solar neigh- 
borhood (Duquennoy and Mayor 1991). The calculations assume a 
temperature of 10 K at 100 AU, as suggested by irradiated disk cal- 
culations (D'Alessio et al. 1998), and a central stellar mass of 0.5Mq. 
The viscosity parameter a has been taken as 0.01, except when stated 
otherwise. The observed values of disk masses and mass accretion rates, 
shown in Fig. 5a and 5b lie within the region bounded by the assumed 
range of initial conditions. Disks with larger initial radii take longer to 
start evolving, since the viscous time scale is tmsc ~ R^ /v oc R. 

The surface density of the similarity solutions behaves with radius 
as 

where Ri and ts are characteristic radius and time (see HCGD for 
details); it goes like oc R~^ at small radii and falls sharply at large 
distances. This last property determines important differences in the 
disk "sizes" measured at different wavelengths, and naturally explains 
the observed disparity between the optical and mm sizes. Disk radii 
measured at millimeter wavelengths are of the order of a few hundred 
AU (Dutrey et al. 1996), while the radii of the disks seen in silhouettes 
in the Orion Nebula cluster at 0.6 /xm are much larger, ^ 500 - 1000 
AU (McCaughrean and O'Dell 1996). 

Figure 6 shows the predicted radii of models in Fig. 5 at 2.7 mm 
and 0.6/im, compared to the observations. Circles indicate the Dutrey 
et al. (1996) observations, while the error bar indicates the range of 
sizes of the Orion silhouettes. To calculate the millimeter sizes the 
two-dimensional brightness distribution of the disk model has been con- 
volved with a Gaussian with the appropriate beam size. Radii at other 
wavelengths correspond to the radii where the optical depth is ~ 1 
at the given wavelength (HCGD). The theoretical predictions compare 
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well with the observations. Since dust opacity increases rapidly towards 
the optical, the outer tenuous regions of the disk can effectively absorb 
background light and produce a large apparent size. In contrast, in 
the millimeter range these outer cooler, low density regions contribute 
little to the surface brightness and the observed sizes are consequently 
smaller (HCGD). The larger millimeter size of the disk in the binary 
(open circle) may be indicative of a circumbinary disk, for which the 
present calculations do not apply. 

Figure 6 shows the predicted size of one of the models at 1.87/im. 
At an age of ~ 0.5Myr, typical of the Orion Nebula cluster, the infrared 
sizes are ~ 20% smaller than the optical sizes, in good agreement with 
observations (McCaughrean et al. 1998). The sharp decline of surface 
density with radius predicted by viscously evolving models naturally 
explains the observed sizes, without the need for heavily truncated 
edges. 

Figure 6 also shows with long dashes the predictions for the ob- 
served sizes in the millimeter range of a disk model with a — 0.001. 
Since tmsc ^ R^/v o: 1/a, disks with small a take much longer to start 
evolving and growing, resulting in sizes much smaller than observed 
at the typical ages of the young population. Thus, measurements of 
disk sizes as a function of time will place important constraints on the 
characteristic value of a in CTTS disks. 

Disk evolution could be considerably different if angular momen- 
tum is lost from outer disk regions through a wind (e.g., Pudritz and 
Norman 1983, Konigl 1989). We have argued elsewhere (Hartmann 
1995) that this is not the case. Similarly, coagulation of disk material 
into bodies that sweep clear the gas will significantly modify this simple 
picture of disk evolution. Nevertheless, it is encouraging that the ob- 
servations can be explained with a viscosity {a ~ 10~^) comparable to 
that estimated in simulations of the Balbus-Hawley magnetorotational 
instability (Stone ct al. 1996; Brandenburg et al. 1996). 

E. Effects of companion stars 

There is a large spread in the mass accretion rates and disk masses 
as a function of age. Some of this range is probably due to uncertain- 
ties in age determinations, errors in accretion rates (for example, due 
to ignoring inclination effects), time- variability, and a range in initial 
conditions. However, the potential effects of companion stars cannot 
be ignored; since at least ~ 2/3 of all systems are binaries (Duquennoy 
and Mayor 1991), it is important to consider the effects of a binary 
companion on the evolution of disk accretion when comparing with 
observations. 

A companion star may prevent the formation of a disk in its im- 
mediate vicinity, and will try to open gaps on either side of its orbit 
(cf. Artymowicz and Lubow 1994, AL; Lubow and Artymowicz, this 
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volume). While the "initial" effects of a binary companion can strongly 
limit disk structure and accretion, it is important to realize that there 
may be secondary effects as well. Specifically, even with a relatively 
distant companion, the inner structure of a viscously-evolving disk will 
eventually be affected by the companion, even though the tidal forces 
are negligible in this region. The reason is that the isolated viscously- 
evolving disk can only accrete if its outer regions expand to take up 
the necessary angular momentum. In the similarity solution described 
above, the mass accretion rate decreases as a power-law in time, because 
as the disk empties out it also expands and thus has an increasingly 
long viscous time. In contrast, if a binary companion limits the expan- 
sion of the disk, once the disk reaches its maximum size, the viscous 
time remains constant, and so the (inner) regions empty out exponen- 
tially with time. (Note that these considerations are relevant only to 
circumstellar disks, not circumbinary disks, which can expand.) 

Figures 5c and 5d show a very simple calculation of this type of 
effect, using the same power-law viscosity used in the standard model, 
but now not allowing the disk to expand beyond a certain outer radius, 
using the boundary conditions discussed by Pringle (1991). One ob- 
serves that when the disk expands to the limiting radius, the accretion 
rate first increases slightly, and then drops precipitously as the disk 
empties out rapidly. 

The significance of this estimate can be seen by noting that the me- 
dian binary separation is roughly 30 AU (Duquennoy and Mayor 1991). 
With the reference disk model used, and estimating a truncation ra- 
dius ~ 1/3 of the binary separation (AL), this would mean that even 
if all binaries originally had circumstellar disks, half of those binaries 
would have their disks empty out by an age of 1 Myr. These estimates 
are roughly consistent with the percentage of ^ 1 Myr WTTS in Tau- 
rus, r^ 45 % (Kenyon and Hartmann 1995); the predicted fraction of 
WTTS could be even higher if indeed the fraction of binaries in Taurus 
is higher than in the solar neighborhood (Simon et al. 1995). 

These estimates of binary effects on disk evolution are rough, and 
the model ignores the effects of the (significant) eccentricities of binary 
orbits (AL), but they serve to illustrate the importance of identifying 
stellar companions to understand disk evolution in individual systems. 
Figures 5c and d do not show a very strong correlation of mass accretion 
rates with binarity, though there is a significant effect on disk masses 
(cf. Jensen et al. 1994; Mathieu et al. 1995; Osterloh and Beckwith 
1995). In any case, many of the systems shown have not been studied 
carefully for potential companions, and in general much work remains 
to be done in this area. 
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IV. SUMMARY AND IMPLICATIONS 



Figure 7 summarizes the ideas presented so far in a sketch of disk 
evolution with time for a single star. After maybe a short initial period 
of high mass accretion rate, the disk remains most of the time in a 
quiescent state. Episodes of high mass accretion rate are triggered 
mostly during the phase where the disk is still immersed in the infalling 
envelope, in which we expect most of the star to be built. After the 
infall ceases and the star emerges as a T Tauri star, the disk evolves 
viscously, M slowly decreases with time, at the same time that the disk 
expands and its mass decreases. 

There are several implications of these results for star and planet 
formation. First, disk accretion during the protostar phase appears to 
be highly variable. This may call into question theories of the birth- 
line, or the initial position of stars in the HR diagram (Stahler 1988; 
Hartmann et al. 1997), which assume steady accretion at the rates of 
infall of the protostellar envelope. It is conceivable that planetesimals 
or other bodies form in the disk during this phase but are swept into 
the star as the disk accretes, perhaps partly accounting for some of the 
accretion variability - our observations really only probe energy release 
in the inner disk, near the star. Second, there appears to be a wide 
variety of disk masses and accretion rates (say, a range of an order of 
magnitude) during the T Tauri phase, produced in part by differences 
in initial angular momenta. This may mean that any consequent plan- 
etary systems which form could have quite different properties. Third, 
the ability of viscous accretion disk models to explain the observations 
so far suggests that substantial migration of material occurs during the 
T Tauri phase; this migration as disks actively accrete may be impor- 
tant in explaining some of the extrasolar planets which lie close to their 
star. Fourth, the presence of binary companions obviously does not al- 
ways prevent disk formation, but they may accelerate (circumstellar) 
disk accretion. Improved mm and submm interferometry, as well as 
infrared speckle searches for companion stars and improved radial ve- 
locity studies to search for close, low-mass stellar companions, will lead 
to a greatly improved understanding of disk evolution. 
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Figure 1. Relationship between accretion luminosity and (a) the ex- 
cess luminosity in the U band, and (b) the luminosity in Br7 for a 
sample of CTTS in Taurus. Data from GuUbring et al. (1997) and 
MuzeroUe et al. (1998c). 
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Figure 2. Histogram showing the distribution of mass accretion rates, 
(a) For T Tauri stars in Taurus and Cha I, with M determined from 
blue spectra or U magnitudes, and (b) For Class 1 sources in Taurus 
and p Oph, determined from Br7 measurements. 
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Figure 3. Relationship between accretion luminosity and bolometric 
luminosity for Class I sources in Taurus (filled circles) and in p Oph 
(open circles) (from MuzeroUe et al. 1998c). 
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Figure 4. Observed mass accretion rate vs. age for CTTS in Taurus, 
Cha I, and p Oph. Mass accretion rates have been obtained by the 
methods described in section IL Ages for the CTTS have been esti- 
mated from the position in the HR diagram and comparison with the 
theoretical tracks from D'Antona and Mazitelli (1994, CMA case). Lu- 
minosities and spectral types were taken from Kenyon and Hartmann 
(1995) for Taurus, Gauvin and Strom (1992) for Cha I, and Greene and 
Meyer (1995) for p Oph. The mean and dispersion of the estimated 
mass accretion rates for Class I sources is also shown for comparison 
(see ILC). The mean age for the Class I sources is assumed to be 0.1 
Myr. 
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Figure 5. Similarity solution for disk evolution, with ly (x R, compared 
to observations. Upper panels: Evolution for isolated disks, (a) Disk 
mass vs. time, (b) M vs. time. Models shown have initial disk mass 
of 0.1 Mq and initial disk radii (marked): 1, 10, and 100 AU. Lower 
panels: Evolution with binary companions, (c) Disk mass vs. time, 
(d) M vs. time. Models are shown for initial disk mass and radius of 
0.1 Mq and 10 AU, and for three binary separations: 30 AU, 100 AU, 
and 300 AU (corresponding to truncation radii of 12, 40, and 120 AU, 
marked). The corresponding evolution for the isolated disk is shown 
for comparison (heavy line). Binaries indicated by open circles. 
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Figure 6. Characteristic disk sizes for viscous evolution as observed 
at 0.6 /im, 1.87 fiia, and 2.7 mm. Models are shown for an initial 
mass of 0.1 Mq and initial radii 10 AU (solid) and 100 AU (short 
dashes). A model with initial radii 10 AU but a ~ 0.001 is shown for 
comparison (long dashes). Data from Dutrey et al. (1996) are shown 
as circles, and the error bar indicates the range of sizes measured in 
the disks seen in silhouettes in the Orion Nebula cluster (McCaughrean 
and O'Dell 1996). Binaries are indicated by open circles. See text. 
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Figure 7. Sketch of disk evolution with time, summarizing the ideas 
presented in this chapter. The disk remains most of the time in a qui- 
escent state, punctuated by episodes of high M as long as the envelope 
feeds mass to the disk. When infall ceases, the disk evolves viscously 
and M slowly decreases with time. 



